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Abstract. This is a continuation of our study about shape parameter, based on an approach very 
different from that of [5] and [5]. Here we adopt an error bound of convergence order O(Vdujd) as 
d —> 0, where < cj < 1 is a constant and d denotes essentially the fill-distance. The constant lu is 
much smaller than the one appears in [5] and [B] where the error bound is O(to^) only. Moreover, 
the constant u here only mildly depends on the dimension n. It means that for high-dimensional 
problems the criteria of choosing the shape parameter presented in this paper are much better than 
those of [5] and [6] . The drawback is that the distribution of data points must be slightly controlled. 
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1 Introduction 

We first review some basic material. The radial function we use to construct approximating functions 
is called shifted surface spline defined by 

h(x) := (-l) m (\x\ 2 + c 2 )? \og(\x\ 2 + c 2 )?, \eZ+, m = 1 + ^, c> 0, x e R n , \,neven, (1) 

where |x| is the Euclidean norm of x, log denotes the natural logarithm, and A, c are constants. 
The constant c is called shape parameter which greatly influences the quality of the approximation. 
Unfortunately, its optimal choice is a big problem and has been regarded as a hard question, not 
only in mathematics, but also in engineering. 

As is well known in the field of RBF (radial basis functions), for any scattered set of data points 
(xi, f(xi)), • • • , (xn, J(xn)), there is a unique function 

N 

s(x) :=*^^Cjh(x — Xj) + p(x) (2) 

3=1 
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interpolating these data points, where Ci,---,cjv are constants to be determined and p(x) is a 
polynomial of degree < m — 1. The only requirement for the data points is that xi, ■ ■ ■ , xn should 
be polynomially nondegenerate. 

The choice of c severely influences the upper bound of \ f(x) — s(x)\. To the author's knowledge, 
there are three kinds of error bound for shifted-surface-spline intcrpolation:algebraic type, exponen- 
tial type, and improved exponential type. Among them the algebraic type shows nothing about the 
effect of c. The exponential type works well for this purpose, as can be seen in j5] and [6] . However 
the improved exponential type works better as will be seen in this article. 

1.1 Function Spaces 

We put restrictions on the approximated functions. 

Definition 1.1 For any a > 0, the class of band-limited functions f in L 2 (R n ) is defined by 

B a := {/ e L 2 {R n ) : /(£) = */ |£| > «t}, 

where f denotes the Fourier transform of f. 
A larger function space is defined as follows. 
Definition 1.2 For any a > 0, 

E a := {/ G L 2 (R n ) : J \f(0\ 2 e^d£ < oo}, 

where f denotes the Fourier transform of f. For each f € E a , its norm is 

i 

\\f\U :={/ |/(£)| 2 e^} 2 . 

Although we only deal with functions from B a or E a in this article, another function space should 
be mentioned. It's denoted by Ch, m and is the so-called native space induced by h. We omit its 
complicated definition and characterization. For these, we refer the reader to [3J [4J [8l [0l [10] . What 
we need here is just B a C E a C C/ ljjn . The proof of E a C Ch, m can be found in [6]. Moreover, each 
function / G Ch,m has a semi-norm denoted by 

1.2 Distribution of Data Points 

The distribution of data points plays a crucial role in our approach. We first review a basic definition 
of i. 

Definition 1.3 Let E be an n-dimensional simplex in R n with vertices V\, ■ ■ ■ ,i> n+ i. For any point 
x G E, its barycentric coordinates are the numbers Ai, • • ■ , A n +i satisfying 

n+1 n+1 

x = XiVi, A, = 1, A; > for all i. 

i=l i=l 

The definition of simplex can be found in [2]. 
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Definition 1.4 For any n- dimensional simplex, the evenly spaced points of degree k are the points 
whose barycentric coordinates are of the form 



(fci/fc, fe/A;, ■ • • , fc I1+ i/fc), ki nonnegative integers and k± + • • • , fc„+i = k. 

As pointed out in [7], the number of such points is equal to the dimension of P£, the space of 
n-dimensional polynomials of degree less than or equal to k. In this paper we use N to denote this 
number. Thus N = dimPu ■ 

In our approach interpolation occurs in a simplex and the centers (interpolation points) are evenly 
spaced points of that simplex. Note that the shape of the simplex is very flexible and hence this 
requirement is not very restrictive. We shall see in the next section that for this kind of interpolation 
there is an error bound which is much better than the case of purely scattered data points, making 
the criteria of choosing c much more meaningful. 

2 Improved Exponential- type Error Bound 

The function h{x) in (1) induces a few basic ingredients of the error bound. First, its Fourier 
transform is 

h{0 = l{\,n)\t\- x - n t^{c\ti\) (3) 

2 

where l(X,n) is a constant depending on A, n|71|U], and K u (t) = t"K. u (t), K, v (t) being the modified 
Bessel function of the second kind pQ. Second, each h(x) corresponds to two constants p and A 
defined as follows. 

Definition 2.1 Let h(x) be as in (1). The constants p and Ao are defined as follows. 

(a) Suppose n — A > 3. Let s = |~ "~ A ~ 3 ] . Then 

s , . (2m + 2 + s)(2m + l + s)---(2m + 3) 
P=1+ 2^T3 and A ° = ^ • 

(b) Suppose n - A < 1. Let s = - [ "~^ 3 1 . Then 
p = 1 and Ao 
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(2m + 2) (2m + 1) • • • (2m - s + 3) ' 
(c) Suppose 1 < n — A < 3. Then 

p = 1 and Ao = 1. 



The following core theorem provides the theoretical ground of our criteria of choosing c. We cite it 
directly from 0. with a slight modification. 

Theorem 2.2 Let h be as in (1). For any positive number bo, there exist positive constants 
Sq, C\, C, U), < w < 1, completely determined by h and bo, such that for any n-dimensional 
simlex Qo of diameter bo, any f G Ch.m, and any < 5 < So, there is a number r satisfying the 
property that ^ < r < bo and for any n-dimensional simplex Q of diameter r, Q C Qo, there is an 
interpolating function s(-) as defined in (2) such that 

\f{x)^s{x)\<c l V~5{u)^\\f\\ h (4) 
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for all x in Q, where C is defined by 

C:=max{v, J-j, p > = i 

where p and c appeared in Definition2.1 and (1) respectively. The function s(-) interpolates f at 
Xi, — ,xn which are evenly spaced points of degree k — 1 on Q as defined in Definitionl.4, with 
k= j. Here \\f\\h is the h-norm of f in the native space. 

The numbers Sq, c\ and ui are given by 5q := 3 c(m+i) w ^ ere 171 appeared in (1), 

ci := y/l(\, n)(2ir)i ^/X'qVzC ^ (I&tt)- 1 

where A is as in (1), l(\,n) appeared in (3), a n is the volume of the unit ball in R n , and A was 

i 

defined in Defi,nition2.1, and ui :— (§) 3C • 

Remark: This seemingly complicated theorem is in fact not so difficult to understand. The number 
5 is in spirit the well-known fill-distance. 

Now, it's easily seen in (4) both c\ and uj depend on the shape parameter c. So does H/IU- In 
order to make (4) useful for choosing c, we still have to convert ||/|U into a transparent expression 
of c. We need two lemmas which we cite from [5] and [6], respectively. 

Lemma 2.3 For any a > 0, f € B a implies f € Ch, m o,nd 

i 

ll/IU < C (m,n) ■ (^) 4 • o 1 -^ ■ c 1 ^ ■ • 

, where 



(27r)-Vm! 
Co(m,n) = . 

Vt(A,n) 

Lemma 2.4 For any a > 0, / € E a implies f € Ch. m o,nd 

n \ 1-"-^ f 1 + n + X c\j\ \(\' , 

ll/IU < a (\,n)c * sup <^ |£| * e ^ )\\f\\ 

, wAere a (A,n) = ^^g^W^i. 
Corollary 2.5 If f E B a , (4) can be converted into 



|/(x)- S (x)|<C BC ^VCe^ ^-J V^||/|| L2 (5) 



where Cb '■— o~ a (2ir) "■y/6na n Ao(167r) _1 rn!. 
Corollary 2.6 If f E E a , (4) can be converted into 



\f(x) - s(x)\ < Cec^VCsxw jiei^e^"^} Q) V*||/|| Bct (6) 



w/iere C B := (2tt)-" v /6na„A (167r)- 1 m!. 
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3 Criteria of Choosing c 



Note that in the right-hand side of both (5) and (6), after every thing independent of c is hxed, 
there is a function of c which may be used to choose the optimal c. However, in Theorem2.2 there 
is still an object dependent of c. That is <5o, the upper bound of 5. For a fixed 5 > 0, decreasing 
c may decrease So and the requirement 5 < S may be violated. For any fixed bo > 0, the minimal 
acceptable c is c = cq := 24p(m + 1)5. Also, C = ^ if and only if c > c\ := Ylpbo- Here ci is 
different from the c\ in (4). 

There is a logical problem about c, 6o, £ an d So. According to Theorem2.2, c appears first, then 
bo, and then So and <5. There is a trick to resolve this logical problem. For any bo > 0, we first let 
C = and <5o := 3 g^ +1 ) = 2(m+i) temporarily. Let (5 < So be fixed. Then (5 < So will always be 
satisfied if c e [c ,oo). After the optimal c is obtained, we let C and So be as in Theorem2.2. The 
consequence is that we can only choose an optimal c from [co, oo). This is a drawback of our theory. 
However, since S can be theoretically arbitrarily small, Co can be very close to zero, theoretically. 

Now, on the right side of either (5) or (6), there is a function of c. Let's call it an MN function 
and denote it by MN(c). Thus for / e B a , 

1 

MN(c) = c'^VCe^ (7) 

, and for / £ E a , 

MN{c) = c^VCs^ n \\^e^-^(^ (8) 

. Our goal is to find c e [c ,oo) which minimizes MN(c). 

In the preceding discussion the domain size bo was fixed. For b not fixed, the way of choosing c 
will be different. We deal with them separately. 



3.1 b fixed 

Note that 



8p 



if cq < c < C\ , 



C 



if Ci < c < oo 



3fco 



, where cq and c\ were defined in the beginning of section3. Thus (7) and (8) can be refined as 



MN(c) = 



y/8pc 4 e 



l T 24pA 



if Co < c < ci 



(9) 



!>0 



-c 4 e 2 



(I) 5 * ifci<c< 



oo 



for / G Bo-, and 

MA^(c) 
for f€E a . 



8£c^ S up^ R n {id^e^-^} (I) 5 * 7 if Co < c< Cl 
3^c 4 sup Cefl „||C| 4 e 2 2 „j.(-) 2a ifci<c<oo 



(10) 
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3.1.1 feB a 

For / g B a , we have the following cases, where k 



£ -4- i!Lt 
2 ^ 24p(5 - 



Casel. 



A - n - 1 > 



and 



/c > 



For any 6q > in Theorem2.2 and positive 5 < 2 (m+i) > ^ 



A — n — 1 > and k > 0, the optimal choice of c for c £ [cq, oo) is to let c = cq := 24p(m + 1)<5. 
Reason: In this case MN(c) in (9) is increasing on [co,oo). 
Numerical Examples: 



Graph of the MN function with 5=0.001 



Graph of the MN function with <5=0.001 



MN(r) 
2.17 




(a) a smaller domain 

Figure 1: Here n = 2, A = 4, a = 40 and b = 10 



0.4 0.6 0.8 

(b) a larger domain 



Case2. 



A - n - 1 > 



and 



k < 



For any bo > in Theorem2.2 and positive S < 



bo 



if 



2(m+l) ' 

A — n — 1 > and k < 0, the optimal choice of c for c G [crj, oo) is c* € [en, ci] which minimizes 
MN(c) of (9) on [c ,ci]. 

Reason: In this case MN(c) in (9) is increasing on [ci,oo). 
Numerical Examples: 
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Graph of the MN function with 5=0.02 



Graph of the MN function with 5=0.02 



MN(c) 



MN(r) 




1.5 
1.0 
0.5 



22 24 26 

(a) a smaller domain 



25 30 35 

(b) a larger domain 



Figure 2: Here n = 2, A = 4, <r = 1 and bo = 2. 
Graph of the MN function with 5=0.01 Gra P h of the MN function with 5=0 01 



MN(f) 



MN(c) 




4.x 10"' 



3.x 10-' 



2.X10" 9 



1.x 10"* 



22 24 26 28 

(a) a smaller domain 



25 30 35 

(b) a larger domain 



MN(f) 



Figure 3: Here n = 2, A = 4, er = l and bo = 2. 
Graph of the MN function with 5=0.009 Graph of the MN function with 5=0.009 

MN(r) 



3.X10" 11 




22 24 26 

(a) a smaller domain 



25 30 35 

(b) a larger domain 



40 





Figure 4: Here n = 2, A = 4, a = 1 and bo = 2. 
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MN(c) 


7.x 10" 18 




6.x 10" 18 




5.x 10" 18 




4.x 10" 18 




3.x 10" 18 




2.x 10" 18 




1.x 10" 18 





Figure 5: Here n = 2, A = 4, a = I and bo = 2. 
Graph of the MN function with (5=0.007 Graph of the MN function with 5=0.007 

MN(c) 



1.5 X lO" 16 



1.x 10" 



5.x 10" 17 - 



22 24 26 28 30 

(a) a smaller domain 




25 30 35 

(b) a larger domain 



Figure 6: Here n = 2, A = 4, er = l and bo = 2. 



All these figures provide a very small part of the entire curve only. In fact the curve increases or 
decreases very rapidly whenever c is far from its optimal value. 



A — n — 1 < 



Case3. 

A — n — 
MN(c) in (9) on [ci,oo). 



and 



k < 



For any bo > in Theorem2.2 and positive 5 < 



2(m+l) ' 



if 



1 < and k < 0, the optimal choice of c G [cq, oo) is the value c* € [ci, oo) which minimizes 



Reason: In this case MN(c) in (9) decreases on [co,ci]. 

Remark: In Case3 if 1 + A — n > 0, MN(c) will be increasing on [c%, oo) and c* = c\. 



Numerical Examples: 



8 



Graph of the MN function with 5=0.02 



Graph of the MN function with 5=0.02 




MN(r) 



0.3 



0.2 



0.1 



24 26 
(a) a smaller domain 



2S 



30 



25 30 35 

(b) a larger domain 



Figure 7: Here n = 2, A = 2,er = l and bo = 2. 
Graph of the MN function with (5=0.01 Graph of the MN function with 5=0.01 



MN(t') 



1.5x10-" 



5.x 10-' 2 






MN(c) 


7.X10 
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-10 
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22 24 26 

(a) a smaller domain 



25 30 35 

(b) a larger domain 



Figure 8: Here n = 2, A = 2, er = 1 and b = 2. 
Graph of the MN function with 5=0.009 Graph of the MN function with 5=0.009 



MN(r) 



MN(t') 




6.x 10' 



4.X10" 12 - 



2.x 10' 



22 24 26 

(a) a smaller domain 



25 30 35 

(b) a larger domain 



40 





Figure 9: Here n = 2, A = 2,er = l and bo = 2. 
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Graph of the MN function with <5=0.008 

MN(c) 



Graph of the MN function with 5=0.008 

MN(c) 




22 24 26 28 30 25 30 35 40 

(a) a smaller domain (b) a larger domain 



Figure 10: Here n = 2,A = 2,er = l and b a = 2. 
Graph of the MN function with 5=0.007 Graph of the MN function with (5=0.007 

MN(c) MN(c) 




22 24 26 28 30 25 30 35 40 

(a) a smaller domain (b) a larger domain 



Figure 11: Here n = 2,A = 2,er = l and b a = 2. 



The following case is quite different from the preceding three cases. For fixed a > 0, the number 
5 cannot be arbitrarily small due to the restriction k > 0. However the optimal choice of c depends 
on the domain size bo. 



Case4. 



A — n— 1 < 



and 



k > 



For any bo > in Theorem2.2 and positive 5 < 2 {m+i) > ^ 
A — n — 1 < and k > 0, the optimal choice of c € [cn,oo) is either c* G [co,ci] or c** € [ci,oo), 
depending on MN{c*) < MN(c**) or AfiV(c**) < MN{c*), where c* and c** minimize MN(c) in 
(9) on [co,ci] and [ci,oo), respectively. 

Reason: In this case MN(c) in (9) may not be monotonic on both [cq, c%] and [c%, oo). 
Remark: In Case4 if 1 + A — n > 0, MN(c) will be increasing on [ci, oo) and c** = ci. 
Numerical Examples: 



10 



MN(r) 



Graph of the MN function with Z? = 10 3 



200 400 600 800 

(a) a smaller domain 



Graph of the MN function with b =10 3 



MN(f) 




1000 



500 



11995 12000 12005 

(b) a larger domain 



Figure 12: Here n = 2, A = 2, a = 1.127 and 5 = 0.03. 



Graph of the MN function with b =lO 



Graph of the MN function with b Q =l0 




(a) a smaller domain 



)0 4000 6000 

(b) a larger domain 



Figure 13: Here n = 2, A = 2, a = 1.127 and S = 0.03 
Graph of the MN function with £> = 10 s 

MN(c) 

MN(e) 



150 



500 1000 1500 

(a) a smaller domain 



Graph of the MN function with £> =10 



1 0.90 
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\ 0.80 




\ 0.75 





5000 10000 15000 20000 

(b) a larger domain 



Figure 14: Here n = 2, A = 2, a = 1.121 and 5 = 0.03. 
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MN(c) 



Graph of the MN function with b =10 6 



MN(r) 



Graph of the MN function with b = 10 6 



10 - 



0.710 - 



MN(c) 



500 600 700 

(a) a smaller domain 



4000 6000 8000 

(b) a larger domain 



10000 12000 



Figure 15: Here n = 2, A = 2, a = 1.127 and 5 = 0.03. 



Graph of the MN function with b =10 



Graph of the MN function with b =l0 7 




650 700 750 

(a) a smaller domain 



1000 1500 2000 

(b) a larger domain 



Figure 16: Here n = 2, A = 2, a = 1.127 and 5 = 0.03. 



3.1.2 feE a 

For / E E a , there are only two cases. 



Casel. 



1 + A - n > 



For any bo > in Theorem2.2 and positive S < 2 (m+i) > if 1 + A — n> 0, the 
optimal choice of c € [cn,oo) is c* € [cn,ci] which minimizes MN(c) in (10) on [cn,ei]. 

Reason: In this case MN(c) is increasing on [ci, oo). Hence its minimum value happens in [en, c{\. 

Numerical Examples: 
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Graph of the MN function with 5=0.01 




4 6 8 

(a) a smaller domain 



Graph of the MN function with 5=0.01 

MN(c) 



8xl0 7 






6xl0 7 






4xl0 7 






2xl0 7 







10 15 
(b) a larger domain 



Figure 17: Here n = 2,A = 2,er = l and bo = 1. 



Graph of the MN function with <5=0.009 

MN(c) 



(1.25 




4 6 

(a) a smaller domain 



Graph of the MN function with 5=0.009 

MN(c) 



8x 10* 






6x 10* 






4x 10 6 






2x 10 6 







10 15 
(b) a larger domain 



Figure 18: Here n = 2, A = 2, a = 1 and 6n = 1. 
Graph of the MN function with 5=0.008 Graph of the MN function with 5=0.008 

MN(c) MN(c) 
0.101- 
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400000 






300000 






200000 






100000 


: 1 





4 6 

(a) a smaller domain 



10 15 
(b) a larger domain 



Figure 19: Here n = 2,A = 2,er = l and b a = 1. 
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Graph of the MN function with 5=0.007 



Graph of the MN function with 5=0.007 



MN(c) 



MN(c) 




1.5 X 10 6 



9 10 11 

(a) a smaller domain 



12 14 16 

(b) a larger domain 



Figure 20: Here n = 2,A = 2,er = l and &o = 1- 



Graph of the MN function with 5=0.006 



Graph of the MN function with 5=0.006 




MN(r) 



10.0 10.5 11.0 11.5 

(a) a smaller domain 



12 14 16 

(b) a larger domain 



Figure 21: Here n = 2,A = 2,er = l and b a = 1. 



20 



The following case happens only when n > 4. 



Case2. 



1 + A — n < For any bo > in Theorem2.2 and positive 5 < 2 (m+i) 7 if 1 + ^ — n < 0, the 



optimal choice of c e [c , oo) is either c* € [co, ci] or c** € [ci, oo), depending on MN(c*) < MN(c**) 
or MN(c**) < MN(c*), where c* and c** minimize MN(c) in (10) on [c ,ci] and [ci,oo), respec- 
tively. 

Reason: In this case MN(c) may not be monotonic on both [co, c\] and [ci, oo). 
Remark: We can apply Matlab or Mathematica to find c* and c** . 
Numerical Examples: 
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Graph of the MN function with 5=0.01 



Graph of the MN function with 5=0.01 



MN(c) 




MN(c) 
8x 10 s - 



6x 10 s - 



4x 10 s - 



2x 10 8 - 



5 6 7 

(a) a smaller domain 



10 15 
(b) a larger domain 



Graph of the MN function with 5=0.009 



Figure 22: Here n = 4, A = 2, er = 1 and b a = 1. 

Graph of the MN function with 5=0.009 

MN(c) 




MN(c) 
0.05 - 



0.03 



5 6 7 

(a) a smaller domain 



10 15 
(b) a larger domain 



Figure 23: Here n = 4, A = 2, er = 1 and b a = 1. 

Graph of the MN function with 5=0.008 

Graph of the MN function with 5=0.008 mn<c) 

5xio 6 - 



4x 10" - 
3x 10 6 - 
2x 10* - 
1 x 10 6 - 



(a) a smaller domain 

Figure 24: Here n = 4, A = 2,er = l and bo = 1 



10 15 
(b) a larger domain 
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Graph of the MN function with 5=0.007 




(a) a smaller domain 



Graph of the MN function with 5=0.007 



MN(r) 



10 15 
(b) a larger domain 



MN(c) 
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0.0020 
0.0015 
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Figure 25: Here n = 4, A = 2, er = 1 and bo = 1. 

Graph of the MN function with 5=0.006 

Graph of the MN function with 5=0.006 MN(c) 
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(a) a smaller domain 



800 



10 15 
(b) a larger domain 



Figure 26: Here n = 4, A = 2,er = l and 6n = 1. 



3.2 



not fixed 



For some domains il the number bo in Theorem2.2 can be made arbitrarily large. For example 
Qo C Q = R n or Qq C Q = {(xx, • ■ ■ , x n )\ < a;, < oo for i = 1, • • • , n}. Such domains are called 
dilation-invariant. In this situation one can keep C = 8p' = ^ in Theorem2.2 by increasing bo, while 
S is fixed. It will make the upper bound (4) better because C can become smaller by increasing c. 

We never decrease bo because it will make both (4) and the upper bound 5q of 6 worse. 

For any S > 0, once the optimal choice of c € [crj,oo) is obtained, we let b Q — Then 

C = 8// = 8 (f ) = and <5 = 3C( ^ +1) = 2 (m+i) ■ Since c > c o = 24p(m + 1)6, we have 
5 < 



3( 8p )(m+i) = 3C(m+i) = an< ^ rec l u i remen ts of Theorcm2.2 are satisfied. 



24p(m+l) 3 (_ 

What's noteworthy is that when bo is increased, in order to keep 8 fixed, one has to add more 
data points to the simplex Q in Theorem2.2. 
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Now, for any 8 > 0, the MN functions become 



MN(c) = y/8p-c^^e 



c e [co, oo) 



for / £ B a , and 



MJV(c) = V^p-c^ - sup ||^|^^e— j ^-J , ce[c ,oo) 

for / e E a . 

Based on (11) and (12), we then have the following criteria of choosing c. 



(11) 



(12) 



3.2.1 / e B a 

Let a > 0, 
Theorem2.2. 



In - 

Let it > 0, o > be fixed, and A; := ^ + • Let / ^ ^ be the approximated function in 



Casel. 



A - n - 1 > 



and 



fc > 



If A — n — 1 > and fc > 0, the optimal choice of c € [cq, 00) is 



c = c := 24p(m + 1)<5. 

Reason: In this case MN(c) in (11) is increasing on [co,oo) 

and 



Case2 
is. 



A - n - 1< 



< 



If A — n — 1 < and fc < 0, the larger c is, the better it 



Reason: In this case MN(c) in (11) is decreasing and MN(c) — > as c — > 00. 

Case3. If A — n — 1 and fc are of opposite signs, MN{c) in (11) is not monotonic and the optimal 
c G [cqjOo) is the number c* minimizing MiV(c), which can be found by Matlab or Mathematica. 
Numerical Examples: 



Graph of the MN function with 5=0.01 



Graph of the MN function with 5=0.01 



MN(c) 



0.4330 - 




MN(f) 
0.530 - 



4500 ^smr 5500 6000 

(a) a = 3.37898 




_2500 3000 3500 4000 

(b) a = 3.37908 



Figure 27: Here n = 2 and A = 2. 
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Graph of the MN function with <5=0.009 

MN(c) 



(a) a = 4.25431 



Graph of the MN function with 5=0.009 



MN(f) 




2 3 

(b) a = 3..95431 



Graph of the MN function with <5=0.008 

MN(c) 



Figure 28: Here n = 2 and A = 2. 

Graph of the MN function with 5=0.008 




(a) a = 4.72359 



2 3 4 

(b) a = 4.42359 



MN(c) 



Figure 29: Here n = 2 and A = 2. 

Graph of the MN function with 5=0.007 
Graph of the MN function with 5=0.007 mn(c) 




(a) a = 5.32697 



2 3 4 

(b) a = 5.02697 



Figure 30: Here n = 2 and A = 2. 
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Graph of the MN function with 5=0.006 

MN(c) 



Graph of the MN function with 5=0.006 

MN(c) 



3.0 



(a) a = 6.13146 



2 3 4 

(b) a = 5.83146 



Figure 31: Here n = 2 and A = 2. 



Note that:(a)in figure27 the two cr's differ only by 0.0001. In Figure28-31 the difference is only 
0.3;(b)in the preceding five examples the optimal choice of c is very sensitive to a, but not to S. 



3.2.2 feE a 

For / € E a , MN(c) in (12) is not monotonic. Hence one should use Matlab or Mathematica to find 
c* G [co,oo) which minimizes MN(c). 



Numerical Examples: 



Graph of the MN function with 5=0.01 

MN(c) 



0.B0 








0.25 








0.20 








0.15 








0.10 










4^"-—. 


6 


i 10 



Graph of the MN function with 5=0.01 

MN(c) 



(a) a smaller domain 

Figure 32: Here n — 2. A = 2 and a = 1 



30000 






25000 






20000 






15000 






10000 






5000 






10 15 
(b) a larger domain 
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Graph of the MN function with 5=0.009 

MN(c) 




4 6 8 

(a) a smaller domain 



Graph of the MN function with 5=0.009 

MN(c) 



1.5X 10 s 






1.0 xlO 8 






5.0 X 10 7 


^ J 





10 15 20 

(b) a larger domain 



Figure 33: Here n = 2. A = 2 and a = 1. 



Graph of the MN function with (5=0.008 




4 6 8 10 

(a) a smaller domain 



Graph of the MN function with 5=0.008 

MN(c) 



1.4X 10 6 






1.2X 10 6 






1.0 xlO 6 






800000 






600000 






400000 






200000 







10 15 20 

(b) a larger domain 



Graph of the MN function with 5=0.007 




0.00015 - 



Figure 34: Here n = 2, A = 2 and a = 1. 

Graph of the MN function with 5=0.007 

MN(c) 



o.s 



7 8 9 10 

(a) a smaller domain 



0.4 



0.2 



10 12 14 16 

(b) a larger domain 



Figure 35: Here n = 2, A = 2 and a = 1. 
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Graph of the MN function with 5=0.006 



MN(c) 
0.00008 h 




10 12 14 16 

(a) a smaller domain 



Graph of the MN function with 5=0.006 

MN(c) 
6 - 

5 1 



3 : 
2 1 



15 20 
(b) a larger domain 



Figure 36: Here n = 2, A = 2 and a = 1. 
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